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Abstract 

We show that the two branches of cohective modes discovered recently in under-doped Cuprates 
with huge spectral weight are a necessary consequence of the loop-current state. Such a state 
has been shown in earlier experiments to be consistent with the symmetry of the order parameter 
competing with superconductivity in four families of Cuprates. We also predict a third branch of 
excitations and suggest techniques to discover it. Using parameters to fit the observed modes, we 
show that the direction of the effective moments in the ground state lies in a cone at an angle to 
the c-axis as observed in experiments. 
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It is generally agreed that the most important theoretical problem posed by the revolu- 
tionary discoveries in the Cuprates is an understanding of the normal states above Tc{x), 
which present new paradigms in the organization of matter. Specifically these states are the 
pseudogap state at low dopings and the non-fermi-liquid state at intermediate dopings. A 
novel broken symmetry with a large order parameter, (estimated jsl to be 0.3jj,B/unitcell 
in Hg(1201) with Tc = QIK), has been discovered in four different families of cuprates [l-5| 
below the pseudogap temperature T*{x), with T*{x) — )■ at x — )■ Xc, the quantum critical 
point. This state can therefore be regarded as a universal feature of the phase diagram of 
the Cuprates. The symmetry of the state with one important discrepancy is consistent with 
the spatial and time-reversal symmetry of a state predicted to arise by formation of ordered 



3] . The marginal non- fermi- 



iquid 



]J] and 



loop currents in a microscopic model for the Cuprates 

properties have been derived from the quantum-critical fluctuations of this state 
the coupling of the fermions to such fluctuations has been shown 15| to promote d-wave 
superconductivity with the right scale of Tc and the superconducting gap. It is therefore 
necessary for a general acceptance of the loop ordered state as the Rosetta stone whose 
deciphering leads to an understanding of all the other universal properties of Cuprates, that 
every aspect of it be convincingly understood. Two crucial aspects towards this goal are 
achieved in this paper. 

If there is an unusual order parameter in the pseudogap state, there should also be 
unusual collective modes there. Recently, inelastic neutron scattering has discovered jsl, ^ 
two branches of weakly dispersive collective modes throughout the Brillouin zone in several 
underdoped Hgl201 for T < T*{x). The integrated strength of these modes at low T is 
estimated to be over 20 times larger than the intensity in the so-called (vr, vr) resonance 8|, |9|. 
Nearly dispersion free modes have already been indirectly inferred through their dominant 
coupling to the fermions in ARPES [l^ and optical conductivity experiments 11 1. The 
energy and the spectral strength of the newly discovered modes make them the most likely 
candidates to explain the ARPES and optical conductivity experiments. 

In this paper we show that three branches of excitations must be present in the loop- 
current ordered state but one of them can not be discovered by neutron scattering. We calcu- 
late the frequency dependence of these modes through introducing the quant um-mechancial 
version of the Ashkin- Teller model which describes the statistical mechanics of the loop or- 
dered state. We predict that the third branch can be discovered through Raman or inelastic 
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x-ray scattering. Moreover, we resolve the major discrepancy of the order parameter found 
in polarized neutron Bragg scattering from the predictions of the classical model. We show 
that the quantum-model gives moments lying on a cone around the c-axis and calculate the 
angle of the cone from the parameters needed to fit the observed collective modes. This 
angle is in agreement with the experiments. We also make several predictions for future 
experiments. 

The four possible " classical" domains of the loop current ordered state are shown in Fig. 
([1]). The "classical" order parameter 12| may be characterized by the anapole vector [l^ L 



/ rfV(M(r) X f) ^ 

J cell 



fi X T/i 



(1) 



where the moment distribution M(r) is formed due to the currents on the four 0-Cu-O 
triangles per unit-cell as shown. An approximate representation of L is given by M(r) ^ 
M^(5^(r — r^), where r^, /i = 1, ..,4 are the location of the four "sites" in any cell at the 
centroid of the moment distribution. These sites are labelled S1,..,S4 in Fig.([T^). The 
orbital moments are either up or down or zero. The four classical domains of Fig.(IT^) 
may be represented by the four values of the angle 6 = 7r/4, 37r/4, 57r/4, 77r/4 that L makes 
with the X— axis. 

The classical statistical mechanics of the Loop-Current state may be derived from the 
Ashkin- Teller model, which is given in terms of a pair of Ising spin per unit-cell ais, Ti^, whose 
eigenvalues, ±1 specify the x and y components of the direction of the vector L and whose 
eigenvectors are |10, 10), |01, 10), |01,01), |10,01); the first two numbers specify eigenstates 
of (Tj3 and the second two of Ti^. 



Hat = ^ ^(o-isCTjs + ^isTjs) + -^CTisTisajsTjs. 



(2) 
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FIG. 1: The four Possible "classical" domains of the loop ordered state are shown. Fig. (e) presents 
the notation discussed in the text for the location of the moments within the unit-cell. 
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Quantum fluctuations among the four possible directions of order together with dissipa- 



tion lead 



14| to a scale invariant spectrum which leads to the observed Marginal Fermi-liquid 



properties in the quantum-critical regime. These quantum fluctuations are specified in the 
basis of the eigenstates of the operator L. The 4 classical loop current states are eigenvalues 
of operators Lj = defined at site i. According to Eq. (8) of Ref.(1l5l). 



[L^ + iLy)\9 >=exp'^\9 > 



The quantum fluctuations are generated from a term in the microscopic hamiltonian: 



(3) 



(4) 



where L^j is the generator of rotations in the space of the four states of Lj, i.e 

U{±^)\Q >= exp{±t^L,i)\e >= \e ± 7r/2 > 
It can then be shown that 



= (l/2)[-(l + 2)t/ + t/' + (-l + z)f/3] 
Ll = 5-2V2(f/ + f/t), 



(5) 



(6) 
(7) 



(8) 



From the operations of f/ -|- f/^ on the basis set chosen in Eq.([2]), it follows that 

/ 1 + i 1 - A 

1 l-i 1+i 

V 2 l-i 1 + i 

\l+i l-i / 

Note that {U + U'') rotates the states \6 > both clockwise and anti-clockwise by tt/2 but 
with matrix elements which are complex. This is represented in terms of the AT model by 
an operator (di -|- cr2 + + ^2)- It is worth noting that in the present problem the kinetic 
energy does not break any lattice symmetry since it is derived from a term L^,- in the original 
model. This is unlike the case of the transverse field Ising model. 

A more convenient representation for calculations on the problem is S't/(4); we introduce 
the direct products, 



(9) 
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The quantum AT model in this representation is 

Hq^t = J2 tiSl + Tl + Sf + T^) + t'{Kl' + KD - J2 ^iSfS^ + TfT^) + ^^f (10) 

i ihj) 

The second term above expresses the effect of f/^ which rotates L by vr. The magnitude of 
At/ J (in the present notation) was estimated [isll] to be of 0(1). It is hard to estimate t', 
which we leave as a free parameter. 

We can parameterize the most general quantum state by the three pairs of angles 6i,(j)i,i = 
1,2,3: 

9/9 9 \ 

1^) = cosy(^cosy|10, 10) + sin-|e^'^^|01,10)j (11) 

+ sin^e^'^i(^cos^|10,01) + sin^e*^3|01,01)) 

Special Case J4 = 0: To understand the results for the general case, it is useful first to 
consider the simple case, J4 = 0, when we simply have two decoupled transverse field Ising 
models. The linear collective modes for the transverse field Ising model were calculated by 
de Gennes [l^. In this case we have (pi = for i = 1,2,3 and ^2 = ^3- Minimizing the 
Hamiltonian we find sin^i = sin ^2 = sin^^ = — t/J. The ground state is 

\G) = (cos^|10) + sin^|01))^(g)(cos^|10) + sin^|01))^ (12) 

There are three eigenvalues at any k, two of which are degenerate, 

ui,2 = 2(J2 - tVk)'/'; /k = (cos(A;,a) + cos(A;,a))/2, (13) 
UJ3 = 4J 

(^1,2 correspond to the propagation of the flipped states in the a and r sectors respectively 
due to the (quantum) transverse-fields, is the energy for flipping the angle of L by vr. 
The modes in this approximation were derived earlier 6| by a simple procedure and called 
the "mixing modes" and the "phase mode" in a semi-classical calculation with no four-spin 
coupling. 

General case: Since the 1 and the 2 components of all operators appear symmetrically, 
(pi = IT /A for all i in the ground state. This leaves 3 real coefficients in the general wave- 
function. Also we are interested only in the case Ji = J2 and the Hamiltonian is invariant 
under interchange of a and r. So the general wave-function, as we have verified by explicit 



calculation, has only two other parameters. It turns out through minimizing numerically 
the ground state energy that the wave-function f|T2l) with just one parameter (beside 0) is 
a good approximate (correct to within a few percent) ground state wave-function for the 
interesting range of parameters. 

To obtain the collective modes in the general case, it is convenient to introduce the 
Holstein-Primakoff transformation generalized to SU{4). We introduce three boson opera- 
tors b,c, and d such that 



= Uc + c^U + + d^b, 
= Ub + b^U + + d^c, 
K^^ = Ud + d^U + b^c + 0^6, 
K^^ = Uc + c^U - bU - d^b, 



= 1- iSc - 2dU 

= 1 - 2b^b - 2Sd 
K^^ = 1- 2b^b - 2c^c 
K^^ = Ub + b^U - c^d - Sc 



(14) 
(15) 
(16) 
(17) 



with U = {1 — b'^b — c^c — d'^d)^^'^. Here b'^ , , d'^ are the creation operators of states |01, 10), 
|10,01) and |01,01) respectively. 

We transform the Hamiltonian to the Boson representation, perform a generalized Bo- 
goliubov transformation and calculate the ground state wavefunction and the eigenvalues 
and the eigenvectors of the excitations. The general results can only be obtained through 
numerical diagonalization of a matrix. We compare such results below with the experiments. 
But first we give the perturbative corrections for J4/J ^ = to the eigenvalues given 
in Eq. (fT3l) for the case J4 = 0: 



4J2 
8J2 
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Note that the degeneracy of the modes uji,uj2 is lifted by J4 7^ 0. This reflects that oper- 
ator K33 couples the a and r modes. For t' = 0, the two modes represent the two linear 
combinations a ± r. For t' 7^ 0, these modes themselves couple and a more complicated 
linear combination varying with k results, as we see through our general results obtained 
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numerically. We can also determine the correction perturbatively in t' for J4 = to get 

t' 



5(a;2)(i)=4j2 



l-(2-3/k)TrTT + /l 



2J2 ' ■'^'2JA 

-l + i2 + 5h)—-U 



t' 



J 

t' 
1 



Now we present the results for the eigenvalues of the collective modes calculated nu- 
merically for any given J, J4,t,t'. For both J4, J > 0, the dispersion of the upper mode is 
smaller than the lower. For J > 0, J4 < 0,t' = 0, the energy of the upper mode at A; = 
progressively decreases as k increases and it must (from the phase diagram of the classical 
AT model) [3] go to at the zone boundary. With t' ^ 0, there is a level repulsion and 
anti-crossing of the two modes. 

Choosing parameters to reproduce the experiments in if (7I2OI with = Q5K, we show 
the results in Fig. ([2]). In the experiments, one of the collective modes is at 40 ± 5 meV 
another one at 50 ± 5 meV at k = 0. The dispersion across the Brillouin zone in both the 
(11) and the (10) directions for the higher energy mode is 5±5 meV, while the lower mode is 
even less dispersive. We also have the very important constraint from the thermodynamics 



of the AT model 



1^ 



191 that to have a transition with no divergence in the specific heat 
— 1 < J4/J < 0, J > 0. In Fig. (E]), the following parameters are used: 

t = 5.5 meV, J = 33.2 meV, J4 = -0.3 J, t' = -1.7t, sinO ^ -0.37 (19) 

We calculate sm6 ^ —0.37 for these parameters. In this case the dispersion width of the 
first mode is about 10% of the energy gap and second one is around 5% of the energy gap. 
The result is shown in figure El These dispersions change with doping and we use a doping 
where a lot of data is available. We can estimate the exact transition temperature Tc from 



using Baxter's [l8| result, exp(— J4/2Tc) = sinh(J/2Tc) to be about 200 Kelvin compared to 
the experimental value of about 250 Kelvin. 

The high energy mode whose eigenvalue is cannot be observed by polarized neutron 
scattering since, it involves change of angular momentum by 2 but should be observable by 
other tools. 

Semi-Classical Angle of the Effective Moments 

The effective direction of the moments is deduced in polarized neutron diffraction experi- 
ments by measuring the intensity of the spin-flipped scattered neutrons at a a Bragg vector 
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FIG. 2: The collective mode dispersions of Ashkin- Teller model with negative as a function of kx 
along the = ky direction. The black curve stands for wi and red curve stands for The 
parameters are taken to fit the experimental results in Hgl201,Tc = 65K.. The blue curve gives 
the calculated third mode. 

Q by sending in the neutron beam at various different polarization directions and using 
standard selection rules [l^]. By this method the direction of M has been determined to be 
45 ± 20 degrees to the normal to the plane both in a sample of Y{123) and if (7(1201). 

The scattering of neutrons from ordered orbital moments leads to a few considerations 
not usually encountered in scattering due to ordered spins. The moments at the four "sites" 
in Fig. (Il^) are orbital moments and so are represented by the transformation properties 
of S* = 1 spins, rather than S = 1/2 spins. The neutron spins are represented by the usual 
spin-1/2 Pauli matrices an, which change the angular momentum of the system by ±1 in 
the spin-flip process. So at any " site" the spin-flp due to scattering of a neutron makes the 
transitions (1 — 0), (—1 — ?■ 0) and vice-versahnt there is no transition from the (1 — > —1) or 
vice- versa. This should be contrasted with ordered spin-systems, for example for S = 1/2, 
where the neutrons scatter from states of 5*^ = — 1/2 to 5^ = 1/2 and vice- versa. 

To calculate the neutron scattering intensity, we must flnd the expectation value of 
(M(Q)) in the ground state. M(Q) is the Fourier transform of the real space magnetic 
moment operator. The representation of M(r) in terms of the order parameter L is ob- 
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tained by inverting Eq. ([T]) to get 

M,(r^) = L, X r^/rl (20) 
Transforming to momentum space 

(M(Q)) = J2 e^^-""' E ^^^^e^«-^ (21) 

where the expectation value is taken in the ground state wave-function, Eq. (fT2|) . which we 
have found to be accurate to a few percent. Using Eq. ([6]) for the matrix representation of 
L, we find 

{L,> = (L„) = -i=r„-'cos9 (22) 

To is the magnitude of measured from the center of the celL Using these we find 

(M,(Q)) = sin ^(^ - 1) (e^Q'i - e^^'^ _ g^Q-a + ^ (23) 

v2 2v2 ^ ^ 

(Mj^(Q)) = -^sin^(^-l)(e*Q-''i+e*Q"2_e*Q.r3_e*Q.r4\ 
v2 2v2 ^ ^ 

Here we have taken the magnitude of M to be unity since we wish to compare with experi- 
ments only the ratio of M in different directions. 

For the momentum transfer Q = (0,a*,c*), as in the experiment and for the domain 
chosen (a in Fig. ([1]) 

(M) = 2isin(a*ry)(^-^sin^(^ -2),0,y2cos^) (24) 

In the other domains M^-, My switch. The tilt angle to the normal is 

-sin^(^ -2)■ 



arctan 



(25) 



2 cos 61 

For fitting the spectrum of collective modes as shown in Fig. ([2]), we needed sin^ ^ —0.37. 
This yields ^ 25 degrees. If we increase t by 20%, the upper mode disperses by about 
10 meV, which is also in the experimental range. Then we can get ^ 33 degreea. The 



experimental range 



9[ is estimated to be 40 ± 20degrees. This means that spin-orbit 



scattering 



20| or additions 2l| to the simplest two-dimensional model may not be required 



to understand the dependence of the intensity of the neutron scattering on the polarization 
direction. 

A prediction, beside that of the third mode which may be discovered through Raman or 
inelastic x-ray scattering, is that the frequency of all three branches of the collective modes 
should — )■ for g — )■ as T — T*{x), because the transition is second order. As already 
noted |3|], the observed ARPES spectra in underdoped cuprates 

To summarize, it is important for an acceptable theory of the Cuprates that the universal 
features of the phase diagram and the principal excitations in each part of it be understood 
consistently and quantitatively from a common set of idea. This paper provides a crucial 
contribution towards this goal through showing that the dominant new fluctuations in the 
pseudogap phase are necessary excitations of the loop-ordered phase which has symmetry 
consistent with the uiversal order parameter observed in underdoped cuprates. An addi- 
tional crucial result is that we find the direction of the effective ordered moments to be 
consistent with experiments through the same quantum-model from which the observed 
collective modes are derived. We also provide several new testable predictions. 
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